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Numerical investigation of supersonic turbulent 
boundary layers with high wall temperature 

By Y. Guo 1 AND N. A. Adams 2 


A direct numerical approach has been developed to simulate supersonic turbulent 
boundary layers. The mean flow quantities are obtained by solving the parabolized 
Reynolds-averaged Navier- Stokes equations (globally). Fluctuating quantities are 
computed locally with a temporal direct numerical simulation approach, in which 
nonparallel effects of boundary layers are partially modeled. Preliminary numeri- 
cal results obtained at the free-stream Mach numbers 3, 4.5, and 6 with hot-wall 
conditions are presented. Approximately 5 million grid points are used in all three 
cases. The numerical results indicate that compressibility effects on turbulent ki- 
netic energy, in terms of dilatational dissipation and pressure-dilatation correlation, 
are small. Due to the hot-wall conditions the results show significant low Reynolds 
number effects and large streamwise streaks. Further simulations with a bigger 
computational box or a cold- wall condition are desirable. 


1. Introduction 

Studies of compressible turbulent boundary layer flows are of fundamental in- 
terest. Understanding of such flows, particularly the effects of compressibility on 
turbulence, may help to improve turbulence models used in many industrial appli- 
cations. Over the years, many experimental efforts have been made (cf. Spina et 
aZ., 1994), and a great deal of knowledge has been obtained about these flows. But 
due to the difficulties encountered in high speed flow experiments, many aspects of 
these flows still remain unclear (cf. Lele, 1994 and Spina et aZ., 1994). Recently 
with the rapid increase of computing power, direct numerical simulation (DNS) has 
become one alternative tool to investigate these flows. In this study, we present our 
first effort in this direction. Preliminary results from three case studies at Mach 
numbers 3, 4.5, and 6 are reported. 

The problem in this study is the compressible turbulent boundary layer flow over 
a flat plate. We consider a region after the leading edge, where the flow is free of 
the shock caused by the leading edge and the boundary layer forms over the wall 
and grows in the streamwise direction. Generally, there are two ways to do DNS 
of boundary layer flows: spatial DNS and temporal DNS (cf. Kieiser &: Zang, 1991 
for a review). A spatial DNS approach (SDNS) is the closest numerical realization 
of a turbulent boundary layer experiment, but it demands an enormous amount 
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of computing resources and therefore is prohibitively expensive, especially in the 
case of compressible boundary layer flows. So far its application has been limited 
to producing a few benchmark data (Spalart & WatmufF, 1993). A temporal DNS 
approach (TDNS), where the flow is assumed to be locally parallel and develops in 
time, is highly efficient and requires no inflow and outflow boundary conditions (cf. 
Kleiser & Zang, 1991). It has been applied successfully to the direct simulations 
of turbulent channel flows (cf. Kim et al } 1987, Coleman et al 1993) and has 
provided considerable details and insights about the turbulence in these flows. But 
due to the use of the parallel flow assumption in its formulation, its applications to 
turbulent boundary layers have suffered certain drawbacks. Firstly, the mean flow 
quantities do not satisfy the Navier-Stokes equations and always exhibit undesirable 
temporal growth. Secondly, the nonparallel effects of a turbulent boundary layer 
cannot be considered. To overcome these problems, Spalart & Leonard (1985) 
and Spalart (1986) assumed self- similarity both for the mean velocity and for the 
Reynolds stresses in their studies of equilibrium turbulent boundary layers and sink 
flows. This assumption is reasonable, but is strictly justified only in the case of sink 
flows. In Spalart (1988), a more general approach was developed using a multiple- 
scale analysis to approximate the local effects of the streamwise growth of the flow. 
Satisfactory results have been obtained for boundary layer flows on a flat plate. 

However, a direct application of Spalart ’s multiple-scale analysis (1988) to com- 
pressible flow is rather tedious and difficult, due to many nonlinear terms in the 
compressible Navier-Stokes equations. In order to avoid this problem, a different 
approach has been used in this study. In the following section, we present the de- 
tails of this approach, together with discussions concerning its application to more 
general cases. Brief comments on the numerical implementation and validation of 
this approach are given in §3. In §4, we present preliminary numerical results of 
turbulent boundary layers at M 0 © = 3, 4.5, and 6. Concluding remarks are given in 
§5. 


2. Numerical method 

It is known that there exist two different scales in a turbulent boundary layer: 
mean flow quantities exhibit large scales and slow variations in the streamwise direc- 
tion, while turbulence quantities are generally of small scale and feist variation (cf. 
Spalart, 1988). In some sense, mean flow quantities can be viewed as global quan- 
tities, carrying upstream information downstream. For example, the downstream 
profiles of mean flow quantities depend on the upstream profiles of these quantities. 
Therefore these quantities cannot be solved by any local procedure without resort- 
ing to any similarity assumption (note that a TDNS approach is a local procedure). 
Turbulence quantities, on the other hand, are generally “local”, being generated 
and dissipated within short time and small space. Intuitively, one would expect 
that it would be more efficient to compute these two different scales with different 
numerical approaches. 

Motivated by the above observation, we developed an approach which computes 
the turbulence quantities locally with a TDNS approach, and computes the mean 
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flow quantities globally by solving the Reynolds-averaged Navier-Stokes equations. 
This approach is formulated in such a way that it allows the solutions of the 
Reynolds-averaged Navier-Stokes equations be marched downstream spatially using 
the Reynolds stresses, which are computed locally at each spatial step by a TDNS 
approach. Special care has been taken to include the nonparallel effects of the tur- 
bulent boundary layer flows into the TDNS approach, which is a natural extension 
of our earlier work on TDNS of growing transitional boundary layers (Guo et a/., 
1994). 

2.1 Solving Reynolds- averaged equations for mean flow quantities within TDNS 

For simplicity, we first demonstrate our numerical approach using an incompress- 
ible boundary layer flow as an example. Application of this approach to compressible 
flow is presented in §2.3. In this report, we designate x as the stream wise direction, 
y as the spanwise direction, and z as the wall-normal direction. Vector notations 
are preferred due to irregular terms introduced by parabolization procedures. 

For any turbulent boundary layer flow, the flow field can be decomposed as 

u = u° + u' , (1) 


where u° represents the stationary part, which usually includes the mean flow, Uq 
(i.e. (0,0) mode), and all stationary streamwise vortices, u$ (i.e. (0,^) modes, 
k y ^ 0). The term u' represents the remaining part which fluctuates with time. 
The governing equations for u° are the Reynolds averaged Navier-Stokes equations 
(in dimensionless form): 

V • (u°u°) = -Vj>° + -^V 2 u° - V • (IPu 7 ) , (2) 

V • u° = 0 , (3) 

where u , u / are the well-known Reynolds stresses, and an overline represents a time- 
averaged quantity. In boundary layers, u° usually has slow variation in the stream- 
wise direction, thus the governing Eqs. (2) and (3) can be parabolized (cf. Fletcher, 
1988 and Guo & Finlay, 1994): the term V 2 is replaced by Vp, 


_2 d 2 d 2 
Vp “ dy 2 + dz 2 ’ 

and the pressure p° in the x — momentum equation is replaced by the mean pressure 
averaged across the y— z plane, p {J. The solutions of the resulting equations then can 
be marched downstream efficiently using a spatial marching scheme. Parabolization 
procedures for more general cases can be found in Fletcher (1988). 

One important issue in solving the parabolized version of (2) and (3) is the numer- 
ical treatment of the Reynolds stress terms. In a turbulent boundary layer, these 
are the dominant terms, and it is desirable to treat these terms implicitly in spatial 
marching schemes in order to avoid any severe restriction on spatial marching step 
sizes. Since the Reynolds stress terms are computed by a TDNS approach in this 
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study, they cannot be expressed explicitly in terms of u°. Thus an efficient way has 
to be found to deal with these terms implicitly. 

It can be shown that in a conventional TDNS approach, the actual equations 
solved for u° are: 

+ Vp • (u°u°) = — Vpp° + -^VpU 0 - Vp ■ (u'u') , (5) 

^ + V P • u° = 0 . (6) 

Note that in (5-6), u° is a function of time and an overline represents a spatial- 
averaged quantity. Compared to the parabolized (2-3), the following term is missing 
from the right hand sides of (5-6): 


Z 0 


(* 01 , * 02 , * 03 , 204 ^ 


( 7 ) 


where 


*01 

*02 

*03 

*04 


du°u° dp$ du'u* 

dx dx dx ’ 

du°v° du'v ' 

dx dx ’ 
du°w° du'w f 
dx dx ’ 

du° 
dx 


Assuming spatial-averaged quantities are good approximations to time-averaged 
quantities, Eqs. (5-6), with the missing term (7) added, constitute a transient prob- 
lem of the parabolized (2-3). The solutions of these equations converge to those of 
the parabolized (2-3) when du° /dt — > 0. Note that the term Zo only acts on the 
governing equations for u°. The streamwise variation of the Reynolds stress terms 
are also included in (7). 

To further demonstrate how Eqs. (5-6) are solved in a TDNS together with the 
term Zo, we discretize the x— momentum equation in the x direction with an Euler 
implicit scheme: 


du°n+l _ d(u°U°) n +I _ d(u°H>°)n+l 1 ( d^n+l 5V«+l\ 

dt dy dz Re \ dy 2 dz 2 ) 

_ a _ (tz°ix°) n+1 - («° u °) n 

\ dy dz ) Ai n+ i 

_ P0n+1 pQn _ (u'tiQn+1 ~ («V)n 
Aln+1 Ax n +i 


where a subscript “n+ 1” denotes the solution evaluated at the downstream location 
x = x„+i, and Ax„+i = x n+ i — x n . At each station x = x n +i, the solution of (8) 
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is iterated in time until a steady-state solution is reached. Then the solution of 
(8) is marched spatially to the next station x = In practice, we find high 

order backward-differentiation-formulas (BDF) with variable step sizes (cf. Hairer 
et al., 1987) are rather satisfactory: these schemes with accuracy up to sixth order 
are unconditionally stable, and the implementation only involves the term Z 0 in an 
existing TDNS code. In our code, a second order BDF is used: 


50 

dx 


^n+1 %n 


1 + 2u n 

1 + W n 


( )n+l ~ (1 + W„)( )„ + 




1 + w„ 


On- 


CS) 


where w„ = (x„+i - x„)/(x„ — x„_i). 

2.2 Computing turbulence quantities with TDNS: nonparallel effects 

Turbulence quantities in boundary layers are not homogeneous locally in the x 
direction. When a TDNS approach is used to compute these quantities, the issue of 
nonparallel effects arises. In Guo et al. (1994), we have discussed extensively this 
issue in the case of transitional flows. Here we adopt a similar methodology. The 
basic idea of this approach is as follows: we assume that the nonparallel flow u can 
be decomposed as the perturbation expansion: 

u = u 0 + Uxe + u 2 e 2 4 , (10) 


where u, (i = 0, 1,2, •••) are periodic functions in both x and y, and e is some 
perturbation parameter. Substituting Eq. (10) into the Navier-Stokes equations, 
we can derive one set of equations for each power i of the parameter e. These 
equation sets all have the periodic functions u, as unknowns. In a sense, the original 
nonperiodic system has been decomposed into periodic subsystems, where then 
Fourier schemes can be used. If Eq. (10) is truncated at i = n— 1, the approximation 
error is given by 

8 = u„e" + u„+ie n+1 + ■ • • . (11) 

If 8 is bounded and decays with n, then Fourier schemes can be used to obtain 
the solution of a nonparallel system. In Guo et al. (1994), we have shown that 
the decomposition (10) can be obtained for transitional boundary layer flows, and 
Eq. (11) is bounded and decays with n. 

For a turbulent boundary layer, we can perform the following analysis. Following 
Spalart (1988), u can be decomposed as 

u = u° + u' = u° + A(x,z)B(x,y, 2 ,<) , (12) 

where u° represents mean flow quantities, and A(x,x) is the “amplitude function”, 
which is proportional to the r.m.s. of u. Both u° and A(x, z) have slow variations in 
the x directions. The term B(x, y, z, t) represents the fast fluctuating part. Taylor- 
expanding u° and A(x, z) in the neighborhood of xo, we have 

1 ^ 2 1|0 

u°(x, z) = u°(x 0 , 2 ) + -fc( x - x o) + 2 j "^f( x _ x °f + ■ " ’ ( 13 ) 
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A(*,*) = A(x 0 ,z) + ^(x - x 0 ) + ^^-(x - x 0 ) 2 + •■■ . (14) 

We denote the length scale of the fluctuating part B(a;, y, 2 , t) as /, and the length 
scale of the mean flow quantities as L. Assuming that d/dx can be scaled with L 
and x — xq with /, we have the expansions (13-14) decay asymptotically at a rate 
0(//Z)"/n!. 

In Spalart (1988), B(:r,y, z,t) was assumed to be periodic in the x direction. 
More generally, we can assume that 

B (x,y,z,t) = B 0 + Bj(x - x 0 ) + B 2 (x - x 0 ) 2 + (15) 


where Bj (j = 0, 1, ■ • •) are periodic in both x and y. Generally, we can assume that 
B, _i is one order smaller than B,. Note that this assumption cannot be proved 
theoretically, but it is less restrictive than assuming B to be periodic. Substituting 
(13-15) into (12), one can show that u' can be written in the form of 

u' = tl(, + Uj(x — X 0 ) + u' 2 (x - Xq ) 2 , (16) 

where u' ( j = 0, 1, - ■ -) are periodic in x and and u( is generally one order larger 
than Then the error defined by (11) is bounded and decays with n. 

Due to the complexity, only the first set of equations (i = 0) are considered in 
this study. The derivative of u is approximated by 


<9u du° Ouq . 

& * aT + 17 + u ' 


(17) 


For a turbulent boundary layer, the nonparallel effect contributed by n[ cannot 
be modeled due to the lack of information about u' x . In a transitional flow, this 
term can be partially modeled, and has been shown to be small in Guo et ai (1994). 
In a standard TDNS approach, the term du° jdx is zero, and thus the nonparallel 
effect contributed by this term is not considered. In this study, the nonparallel 
contribution of this term is considered by adding the following term to the right- 
hand side of the governing equations solved by a TDNS approach, 


Zi 


(* 11 , 212 , * 13 , * 14 ) 


T 

, 


(18) 


where 


z n = -2(u-u 0 )^ , 

( tK^° / 0x0“° 


. dw° 


du° 


Z„ = -(U-U) w -(W-W)-^, 
*14 = 0 . 


The term Z\ is equivalent to the forcing term Zi in Guo et ai (1994), and it only 
acts on the governing equations for u'. 
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2.S Application to compressible boundary flow 

In the case of compressible boundary layer flows, the governing equations we have 
solved can be written in the dimensionless conservative form as: 


au 

dt 


dF dG 
dx ^ dy 


+ ^ + Z 0 + Z 1 
dz 


(19) 


The vector U is the solution vector defined by U = (p,pu,pv, pw y E) T , where 
p is the density normalized with the free-stream density p and u, u, and w are 
the velocity components normalized with the free-stream velocity u^. Here, a 
subscript w oo” denotes the quantity evaluated at the freestream and an asterisk 
denotes a dimensional quantity. The total energy E is defined by 


E = -1— -p + ^(u 2 + v 2 + w 2 ) 

K — 1 2 


where k is the ratio of specific heats and p is the pressure normalized with p* 0 qU^ 2 . 
The pressure p is related to the temperature T by a perfect-gas law 

= pT . 


Sutherland’s law is used for the viscosity calculation. The vectors F, G and H are 
the standard flux vectors in the x, ?/, and z directions respectively. 

The term Zo is derived in the same way as described early in §2.1. It has the 
following form: 




—pu 

-pu 2 - pg + (T«)r 

-puv + (r^) T I + Zr, 

~Ww + (r^7)r 

-u(E + p)~ (g7)r + + (vT zy )r + (t5^7)r ■ 


( 20 ) 


where r n , and t xz are the components of the shear stress tensor and q x is the 
heat flux component in the x direction. Here an overline denotes a spatial- averaged 
quantity. A subscript U T” denotes the quantity been evaluated in a TDNS scheme. 
The term Z m contains the terms missing in the terms ( )x> due to the use of Fourier 
schemes. It has a form of 


Zm — (*i m ? ^ 2m j ~3m > ^4 m •> m ) 


( 21 ) 


where 


Z\m — 0 , 
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d / ji \ dw 2 d / p d u\ 2 d f p du\ 

dx \ Re) dx 3 dy \Re dx ) 3 dz \ Re dx ) 

_ d / ft \ dv 
dx\Rt/dx 



It can be shown that Z m is very small, and can be neglected. For completeness, we 
retain this term in our scheme. 

Note that in Eq. (20), a parabolization procedure similar to the one described in 
Fletcher (1988) has been used: the pressure p in the x—momentum equation is first 
replaced by the mean pressure p®; secondly, the mean pressure pj} in the subsonic 
section of the boundary layer is extrapolated from the adjacent supersonic layer. 
For a boundary layer over a flat plate, we set dp^/dx to be zero. 

The term Z\ in Eq. (19) is given by 


Zi 


(*11? *12? *13, *14, *is) 7 


( 22 ) 


where 


- “°>lr _ ipu ~ ' >u); sx ■ 


* 11 = 

*12 = 

*13 = 

*14 = ~(U “ W U ) 
*15 


du° dp° d(p 2 u°) 




dx 


du° 

dx 

o,dpw — du° 

aT 97 


-(« - + p°)-(e+p-e°- p° 


2.4 Some comments 

It can be shown that the governing equations obtained from our approach for 
an incompressible boundary layer over a flat plate axe identical to those of Spalart 
(1988) when the coordinate transformation in Spalart (1988) is not used. In a sense, 
our approach is similar to Spalart 7 s (1988) in that in both two scales are utilized. In 
our approach, the long-scale analysis is explicitly applied to the governing equations 
for mean flow quantities (i.e. (0, k y ) modes), and the short-scale analysis to the 
governing equations for turbulence quantities. When there are large streamwise 
streaks, e.g. in a three-dimensional boundary layer, our approach may be more 
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general since the related governing equations for mean flow quantities are explicitly 
parabolized, making the use of spatial marching schemes more justifiable. 

In Spalart (1988), coordinate tilting has been introduced to model the nonparallel 
effects of boundary layers. In Guo et al (1994), our experience with transitional 
compressible boundary layer flows shows that the coordinate transformation tech- 
nique does not guarantee an improved result, partly due to the difficulty in selecting 
a proper tilting angle. In this study, this technique is not used. 

Compared to the parabolized stability equations (PSE) approach developed for 
transitional flows (Bertolotti, Herbert & Spalart, 1991), we have the following com- 
ments. The PSE has been proven to be a very powerful tool in studying the transi- 
tion process in spatially growing boundary layers up to the early stage of breakdown. 
The possibility of developing a PSE type approach for turbulent boundary layers 
has been appealing (Zang, 1991). To our understanding, the PSE relies on the 
spatial marching of shape functions (eigenfunctions) of all modes in the streamwise 
direction. This requires the spatial evolutions of all shape functions to be slow in 
this direction. In turbulent boundary layers, only the mean flow u* and the sta- 
tionary spanwise modes have slowly evolving shape functions in the streamwise 
direction. The shape functions of other modes do not have slow variations in the 
streamwise direction, and thus cannot be marched downstream by a spatial march- 
ing scheme. The method we described in this study has explored the possibility of 
marching downstream the shape functions of Uj and u§, and computes the other 
modes locally. Thus it can be viewed as a version of PSE for turbulent boundary 
layer flows. 

3. Implementation and code validation 

The numerical method discussed in §2 has been implemented in a TDNS code 
developed by Adams (1993). Fourier collocation methods are employed for the 
spatial derivative calculations in the spanwise and streamwise directions, and a 
compact central finite-difference (Pade) scheme is employed in the non-periodic 
wall-normal direction. The solution is advanced in time fully explicitly by a third- 
order low- storage Runge-Kutta scheme. More details of this TDNS code can be 
found in Adams (1993). 

In order to maintain flexibility, both streamwise and spanwise spatial-averaging 
can be performed on the quantities at x n +i. For the quantities at x n and r n _i, 
time-averaging is also implemented as an option, since the spanwise modes in 
some cases evolve with time t and a time- averaging is needed to obtain stationary 
solutions. In this study, spanwise spatial-averaging is not used. 

In order to validate the code, we first simulate a laminar boundary layer flow at 
Moo = 4.5, with all non-zero streamwise modes set to be zero (i.e. the Reynolds 
stress terms are set to be zero). Similarity solutions for the boundary layer were 
used as the initial profiles to start the downstream spatial marching. At each spatial 
step, the solution is advanced in time until a steady-state solution is reached. Once 
the steady-state solution is achieved, the solution is marched downstream again one 
step. At all spatial steps, our code can duplicate the similarity solutions accurately. 
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Simulations of turbulent boundary layer flows cannot be validated directly due 
to the lack of comprehensive database for these flows. In the following section, 
we present preliminary results for the turbulent boundary layers at three different 
Mach numbers. 

4. Preliminary results 

4.1 Simulation parameters and procedures 

In order to study the compressibility effects at different free-stream Mach number, 
we have done three case studies at Mqq — 3, 4.5, and 6 . The simulation parameters 
for these three cases are listed in Table 1 . The sizes of the computational boxes, L x 
in the x direction, L y in the y direction, and L z in the wall-normal direction, are 
given in the terms of both boundary layer thickness, 60 , and wall unit. The numbers 
of grids in the x, y, and z directions are given by 1ST*, N y , and N z respectively. 
Sampling time t 9 is non-dimensionalized by the free-stream variables. In all cases, 
the laminar adiabatic wall temperatures were used as the constant wall temperatures 
in the simulations, with T ^ = 61.15A\ 

TABLE 1 . Simulation parameters. 


M 00 

Reg 

,/Too 

L z 

Ly 

L z 

N x xN y xN z 

ts 

3 

3015 

2.5 

1.676 0 

527+ 

0.966 o 

300+ 

3.596 0 

1133+ 

192x144x180 

34 

4.5 

2618 

4.38 

1.596o 

260+ 

0.95<5o 

155+ 

3. I 860 
520+ 

180x144x190 

32 

6 

2652 

6.98 

2.166o 

229+ 

1.36o 

137+ 

3.686o 

390+ 

180x128x200 

19 


The computation at M = 4.5 was first performed. The TDNS of a subhar- 
monic transition process at Moo = 4.5 (Adams, 1993) was continued until the flow 
has gone through the breakdown stage and become turbulence. Then the TDNS 
transition data were used as the initial solution to start the spatial marching with 
the approach described in § 2 . Since the TDNS transition data are not the solutions 
of the Reynolds- average Navier-Stokes equations, spatial transients are expected in 
the first few steps. Usually, the transients diminish after a few steps. The spa- 
tial extent affected by the transient solutions depends on both the character of the 
Reynolds- averaged Navier-Stokes equations and the initial conditions. This issue is 
beyond the scope of this report. 

All computations are performed on a NEC-SX3 computer with one processor, 
with the code running typically at 1.8 GFLOPS. To cut down computation cost, a 
smaller box was used at the first step, and the solution is advanced in time until 
a statistically steady-state is reached. Then the solution was marched downstream 
one step spatially with a bigger computation box. In the cases of = 3 and 6 , 
the solutions from the Moo =4.5 case were first rescaled and used as the initial 
data to start the spatial marching. Similar to the case of Moo = 4.5, smaller boxes 
were used in the first two steps, and then a bigger box was used at the third step. 
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FIGURE 1. (a) Time history of skin friction coefficient, C/. Legend: step 1; 

step 2; step 3; step 3 with a big box; o sampling period. ( b ) Mean 

profiles of u, T and p at Moo = 3 over t = 438 - 472. The profiles at t = 472 are 
indicated by . 




FIGURE 2* One-dimensional spectra of u' at z = 30 + and Moo = 3: (a) in the x 
direction; ( 6 ) in the y direction. Legend: u\ v ; w\ T; p. 

Fig. 1(a) shows the time history of skin coefficient, C /, during the first three steps 
at Moo = 6. Here t is non-dimensionalized by the free-stream quantities. It can 
be seen that at each station, Cf is able to reach a statistically steady state. The 
data presented in this report is taken from the third step. At the first step, a large 
variation of C / is observed. We believe this is due to the large transient caused by 
the initial conditions. Note that the initial data is rescaled from that of M 0 0 = 4.5 
case. A similar pattern has also been observed at Moo = 3. 

Fig. 1(6) shows the variations of the mean profiles, u, T and ~p over the time 
period 438 — 472 at Moo = 3. It can be seen that though the boundary layer does 
not grow with t, it does vary with t. We expect that the fluctuations of the mean 
flow quantities with time will become smaller when a bigger computational box is 
used. 
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In all cases, one-dimensional spectrum is monitored to ensure that the computa- 
tions are well resolved. Fig. 2 shows the one- dimensional spectra near the wall at 
Moo — 3. Here the wavenumbers, k x and k y , are normalized by the boundary layer 
thickness. It can be seen that in both plots, the difference between maximum and 
minimum spectral components is maintained to be at least 10 3 5 to 10 4 . Similar 
situations are found in the cases of Moo = 4.5 and 6. 

Figs. 3-5 show the two-point correlations for all three cases. Compared to incom- 
pressible channel or boundary layer simulations (Kim et a/., 1987 and Spalart, 1988), 
the values of these functions are too high near the centers of the boxes. A similar 
situation was also reported in the compressible channel simulations of Coleman et 
ai (1993). In Coleman et ai (1993), acoustic effects are found to be responsible 
for high correlations observed in the compressible channel flows. Here this issue has 
not been investigated thoroughly so far. In all our cases, rather large streaks have 
been observed (cf. §4.2). We speculate that the existence of these large streaks near 
the wall might be the primary reason for the high two-point correlations. 

Note that the downstream developments of the solutions have not been investi- 
gated in this study, since such investigations require the solutions to be marched 
downstream a few steps. So far we have only marched the solution downstream 
three steps at most. The solutions at these steps might still be affected by the 
spatial transients. 

4-2 Large streamwise streaks 

In all cases, large streamwise streaks have been observed. Fig. 6 shows the contour 
surface of u — 0.4 at M 00 = 3. In the figure, the spanwise dimensions of the streaks 
are about 100 + , which are close to those found in an incompressible boundary layer 
flow. The streamwise dimension of these streaks cannot be obtained from the figure 
since the computational box is not large enough. Fig. 7 gives an impression of the 
extent of these streaks in the wall-normal direction. Compared to the boundary 
layer thickness (<5o = 314 + ), they are rather large, affecting almost one third of the 
boundary layer. 

Figs. 8 and 9 show a similar situation at Moo = 4.5, except that here the streaks 
are even bigger compared to the boundary layer thickness (£0 — 163 + ). So far we 
have not been able to study the spanwise and streamwise extents of these streaks 
systematically, primarily due to the their large size together with the present lim- 
itation to small computational boxes (in terms of wall units). Since these streaks 
are rather large compared to the computational boxes, their evolution affect both 
mean flow quantities u° and turbulence quantities u' (note that u' is solved from 
(19) correctly only if u° is steady). Thus the quality of our DNS data needs still 
to be assessed, and readers should be aware of that in comparing our DNS results 
with other sources. 

The mechanism for the formation of these large streaks are not clear to us. In 
compressible channel flows, DNS results have shown that compressibility effects may 
be responsible for rather large structures in the streamwise directions (Coleman et 
a/., 1993, Coleman, 1993 and Moser, private communication, 1994). In our cases, we 
speculate that both compressibility effects and low Reynolds number effects due to 
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FIGURE 3. Two-point correlations at z — 1S.8 + and Moo = 3: (a) in the x 
direction; (6) in the y direction. Legend: u\ v; w\ T; p. 



FIGURE 4. Two-point correlations at z = 17. 3 + and Moo = 4.5: (a) in the x 
direction; (6) in the y direction. Legend is the same as Fig. 3. 
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FIGURE 5. Two-point correlations at z = 15.6 + and M 0 0 = 6: (a) in the x 
direction; (6) in the y direction. Legend is the same as Fig. 3. 
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FIGURE 10. Profiles of (a) mean temperature, T, and ( b ) mean kinematic viscosity, 
v. Legend: Moo = 3, Moo — 4.5, Moo = 6. 



FIGURE 11. Mean streamwise velocity, a, after Van Driest transformation. Leg- 
end: Moo = 3, Moo = 4.5, - - Moo = 6, 2.441n(y+) + 5.5. 

hot-wall conditions may play important roles. With an increase in wall temperature, 
the density of air decreases and its viscosity increases. This effectively reduces the 
Reynolds number in the near wall region. Fig. 10(a) shows the averaged temperature 
profiles across the boundary layers at different Mach numbers. The corresponding 

profiles for the kinematic viscosity, F, are shown in Fig. 10(6). It can be seen 

"—2 

that the increase of v is roughly proportional to T . Therefore, the reductions of 
effective Reynolds numbers in the near wall regions due to the hot-wall conditions 
are significant. 

4.8 Mean flow quantities 

Fig. 11 shows the streamwise velocity profiles after the Van Driest transforms. 
Here the Van Driest transformation is defined by 
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It can be seen that in all three cases, there are no apparent log-law regions. Note that 
at similar Reynolds numbers, some log-law regions can be observed in incompressible 
boundary layer flows (cf. Spalart, 1988). Thus we may conclude that the low 
Reynolds number effects due to the hot- wall conditions in our cases are significant. 

Table 2 shows the skin friction coefficients, in comparison with the results of Van 
Driest transformation II, computed by Zeman (1993) at Re$ ~ 10000. Note that 
the Van Driest transformation II is usually formulated in terms of Re x . Due to the 
difficulty in converting Re x to Re$ for a compressible turbulent boundary layer, we 
did not compute Cf at the relevant Re $ so far. 

TABLE 2. Computed skin friction coefficients compared to the calculations of 
Zeman (1993) using Van Driest transform II at Re$ = 10000. 


Moo 

C/xlO -3 

c f /c fi 

Van Driest II 

3 

1.854 

0.522 

0.58 

4.5 

1.729 

0.470 

0.40 

6 

1.110 

0.303 

0.30 


Due to significant low Reynolds number effects in our DNS data, the effectiveness 
of Van Driest transformation cannot be demonstrated in this report. Readers are 
referred to Huang & Coleman (1994) for a further discussion on this issue. 

4-4 Compressibility effects 

Fig. 12 shows the profiles of turbulent Mach numbers. Note that Alt does not 
increase to the value that one would expect at M & = 6. We believe that this 
is primarily due to the higher wall temperature at Moo — 6. We have seen in 
Fig. 10(6) that the kinematic viscosity, 17, is roughly proportional to T . A lower 
than expected M t at M 0 0 = 6 is presumably caused by the high 17, which damps 
out the fluctuations in the near wall region. 

In homogeneous flows, compressibility effects in the turbulent kinetic e nergy 
(TKE) equation can be measured by the dilatational dissipation, 'ptd = , 2 , 

and the pressure-dilatation correlation, pd = pu\ ^ (cf. Friedrich, 1993, Sarkar,1993, 
Spina et a/., 1994 and Lele, 1994 for a review). Recent work by Huang et al (1994) 
also indicates that this is true for the near-wall turbulence. The inhomogeneous 
terms resulted from the compressibility of fluid are small and confined to the region 
very close to the wall. In Fig. 13, we have plotted these two terms, in comp arison 
to the incompressible part of the dissipation, solenoidal dissipation ~pe a — Jtu /u?-, 
where is the component of vorticity. It can be seen that in all three cases, these 
two terms are small compared to solenoidal dissipation, ~pc 3 . This indicates that the 
influence of fluid compressibility on TKE is rather small in the flat plate bound- 
ary layer flows investigated here. Similar observations are also made by Huang et 
al. (1994) for near- wall turbulence in compressible channel flows at Mach number 
3, and by Dinavahi &z Pruett (1993) and Adams (1993) for late stage transitional 
boundary layer flows at Moo = 4.5. 
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FIGURE 12. Turbulent number Mach number, Mi = M^yJ u”u” j\/T. Legend: 
Moo = 3, Moo — 4.5, Moo = 6. 




FIGURE 13. Compressibility effects: — solenoidal dissipation, pe 3 , di- 

latational dissipation, pQ, pressure-dilatation correlation, pd . (a) M 0 0 = 3. 

'pea and pd have been enlarged 100 times. ( b ) = 4.5. and pd have been 

enlarged 10 times, (c) Moo = 6. ped and pd have been enlarged 10 times. 
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FIGURE 14. Favre- averaged Reynolds stress terms normalized by wall unit: 
u n u n \ t/V'; w n w n ) and u lt w n . 



FIGURE 15. The r.m.s. of turbulence quantities: u* , v ' , u> # , 

r, - - p'. (a) Moo = 3; (6) Moo = 4.5; (c) = 6. 
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Log (ig 

FIGURE 16. Spectra of u in the x direction taken from z/6q = 0.04 
0.1 , 0.16 , 0.42 , and 0.85 — at Moo = 4.5. 



FIGURE 17. Kolmogorov scales, normalized by wall unit. Legend: M 0 0 = 

3, Moo — 4.5, Moo = 6. 

j.5 Other quantities 

Since our DNS data are severely affected by low Reynolds number effects, large 
streaks, and small computational boxes, we do not expect to get good statistics in 
all three cases. In this subsection, we shall only present a few basic quantities to 
give readers a qualitative impression. 

The Favre-averaged Reynolds stress terms normalized by the wall units at Moo = 
3 are shown in Fig. 14. Fig. 15 shows the turbulence fluctuation u' for the three 
cases. It can be seen that the density fluctuations, p f , have larger amplitudes near 
the edges of the boundary layers in all three cases. 
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5. Discussions 

In all our computations, the computational boxes are rather small in terms of 
wall units, but not small at all in terms of boundary layer thickness. That is 
due to the fact that the hot-wall conditions increase the wall units in terms of 
boundary layer thickness. Apparently for the cases studied in this report, larger 
boxes are needed in order to give improved turbulence statistics. This, however, is 
expected to be computationally very expensive if the resolution used in this study 
is to be kept. One may argue that coarser grids might be sufficient since the hot- 
wall conditions reduce the effective Reynolds number in the near wall region and 
turbulence scales are large away from the walls. This argument is not supported by 
our numerical experience so far. Fig. 16 shows the one-dimensional spectra taken 
at different distances from the wall at Moo = 4.5. It can be seen that the spectra 
have almost the same value at high wavenumber from z/6q = 0.21 to 0.85. This 
indicates that the resolution requirement is not reduced at all away from the wall. 
In order to find the explanation, we plot in figure 17 the Kolmogorov scale, /^, as 
functions of z. In the figure, we can see that /j!" has roughly the same value in 
most part of the boundary layer. Similar situations can be seen in Moo = 3 and 6 
cases. Since one needs rather dense grids in the near wall region to resolve mean 
flow profiles, the resolution requirement may be more or less the same throughout 
the whole boundary layer. This observation is basically supported by the one- 
dimensional spectra in Fig. 16. Apparently, this situation is rather different from 
that in incompressible boundary layer flows, where dense grids are needed only in 
the near wall region. From the above discussion, we can conclude that turbulent 
boundary layers with hot-wall conditions require a larger computational effort than 
in the cold-wall case. 

In Fig. 17, one can also see that is much smaller than that in an incompressible 
fluid case, where / jj" % 2 + (Moser, private communication, 1994). This suggests that 
the grid spacing rules derived from DNS of incompressible boundary layer flows 
cannot be applied to a compressible flow case. 

6. Concluding remarks 

In this report, we have presented a numerical approach for DNS of compress- 
ible turbulent boundary layer flows, which solves the parabolized Reynolds- average 
Navier- Stokes equations for mean flow quantities and computes turbulence quanti- 
ties locally with a temporal DNS approach. Preliminary results from three cases 
at Moo “ 3, 4.5, and 6 are presented. All our results indicate that compressibil- 
ity effects on turbulent kinetic energy equations, i.e. dilatational dissipation and 
pressure-dilatation correlation, are small. Due to the hot-wall conditions used in 
all cases, these preliminary results are severely affected by low Reynolds number 
effects, large streaks, and small computational boxes. Hot-wall conditions seem to 
be more difficult to deal with than cold-wall conditions. Grid spacing rules derived 
from the DNS of incompressible boundary flows are not valid in a compressible case. 
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